Temperature dependent resistivity in bilayer graphene due to flexural phonons 
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We have studied electron scattering by out-of-plane (flexural) phonons in doped suspended bilayer 
graphene. We have found the bilayer membrane to follow the qualitative behavior of the monolayer 
cousin. In the bilayer, different electronic structure combine with different electron-phonon coupling 
to give the same parametric dependence in resistivity, and in particular the same temperature T 
behavior. In parallel with the single layer, flexural phonons dominate the phonon contribution to 
resistivity in the absence of strain, where a density independent mobility is obtained. This contri- 
bution is strongly suppressed by tension, and in-plane phonons become the dominant contribution 
in strained samples. Among the quantitative differences an important one has been identified: room 
T mobility in bilayer graphene is substantially higher than in monolayer. The origin of quantitative 
differences has been unveiled. 

PACS numbers: 72.10.-d, 72.10.Di, 72.80.Vp, 81.05.Uw 
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I. INTRODUCTION 

Bilayer graphene continues to attract a great deal 
of attention because of both fascinating fundamental 
physics^ and possible applications.— Recent realization 
of suspended monolayeii^ and bilayer- graphene sam- 
ples made possible a direct probe of the intrinsic, un- 
usual properties of these systems. In particular, intrin- 
sic scattering mechanisms limiting mobility may now be 
unveiledi^ It has been recently shown that in suspended, 
non-strained monolayer graphene room temperature T 
mobility is limited to values observed for samples on 
substrate due to scattering by out of plane — flexural 
— acoustic phononsj^ This limitation can, however, be 
avoided by applying tension. Bilayer graphene has a dif- 
ferent low energy electronic behavior as well as different 
electron-phonon coupling. It is then natural to wonder 
what is the situation in the bilayer regarding electron 
scattering by acoustic phonons, and in particular by flex- 
ural phonons (FPs). 

In the present work the T dependent resistivity due 
to scattering by both acoustic in-plane phonons and FPs 
in doped, suspended bilayer graphene, has been investi- 
gated. We have found the bilayer membrane to follow 
the qualitative behavior of the monolayer parent j^ii Ex- 
plicitly, at experimentally relevant T, the non-strained 
samples show quadratic in T resistivity with logarithmic 
correction, g ^ a-nd constant mobility. Electron 

scattering by two FPs gives the main contribution to the 
resistivity in this case, and is responsible for the de- 
pendence. Suspended samples may also be under strain 
either due to the charging gate^ or due to the experi- 
mental procedure to get suspended samples, or even by 
applying strain in a controlled wayj^^^^ Under uniaxial or 
isotropic strain u the T dependence of resistivity due to 
FPs becomes quartic at high strain u ^ u* , g ^ T'^/v?, 
and quadratic at low strain u ^ u* , g ^ T'^/u, where 
u* ~ 10~'*r(K). These contributions are weaker than 



that coming from scattering by in-plane phonons, and 
in strained samples the latter dominates resistivity, as 
has been found for monolayer graphene;^ An interesting 
quantitative difference with respect to suspended mono- 
layer has been found. In the latter, room T mobility /i is 
limited to values obtained for samples on substrate due 
to FPs, ^ 1 m^/(Vs)i^ In bilayer, quantitative differ- 
ences in electron-phonon coupling and elastic constants 
lead to room T enhanced /i ^ 10 — 20 m^/(Vs), even in 
non-strained samples. 

The paper is organized as follows. In Sec. |TT]we intro- 
duce long wavelength acoustic phonons in the framework 
of elasticity theory. We show how the dispersion rela- 
tion of FPs are affected by the presence of tension over 
the sample. Then, we review the electronic low-energy 
description of bilayer graphene and deduce the electron- 
phonon coupling within this approach in Sec. IIIII The 
variational approach used in order to study the T depen- 
dent resistivity due to scattering by in-plane and FPs 
and a summary of our results in different regimes of T 
and strain are presented in Sec. IIVI Sec. |V] is devoted 
to discuss the implications of these results, the differ- 
ences between monolayer and bilayer and some experi- 
mental consequences. Finally, we expose our conclusions 
in Sec. IVII Some technical aspects are treated in detail 
in appendices. In Appendix. |X] we present the collision 
integral due to scattering by acoustic phonons, and in 
Appendix [B] its linearized form is derived. Details on the 
calculation of the resistivity using the variational method 
are given in Appendix [C] In Appendix |D] we discuss how 
anharmonic effects are partially suppressed by the pres- 
ence of strain. 



II. ACOUSTIC PHONONS 

At long- wavelengths the elastic behavior of monolayer 
graphene is well approximated by that of an isotropic 
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continuum membraneiiii^ whose free-energy reads^i^iii with a = . The typical value is a ~ 4.6 x 10~'^m^/s 



dxdy{\u%+2^iuiA. (1) 



The first and second terms in Eq. ([T]) represent the bend- 
ing and stretching energies, respectively. Summation 
over indices is assumed. In-plane distortions are denoted 
by u(r) and out-of-plane ft.(r), with r — (x, y), such that 
the new position is X{y) = (x, y, 0) 4- [wa;(r), Uj^(r), /i(r)]. 
To lowest order in gradients of the deformations the 
strain tensor appearing in Eq. ([1]) is 



1] — 2 [diUj 



(2) 



and owing to the same argument the factor ^ = 
y/1 + |V/iP in the measure is neglected. The parame- 
ter K is the bending rigidity and A and fj. are in-plane 
elastic constants. Typical parameters for graphene are 
K « 1 eV and ^ ~ 3A « 9eVA"^^^>i^ with mass den- 
sity p = 7.6 X 10~'^kg/m^. 

In the case of bilayer graphene, as long as we are not 
at too high T to excite optical phonons we may, from the 
elastic point of view, regard bilayer graphene as a thick 
membrane with mass density and elastic constants twice 
as high as those for single layer 

The dynamics of the displacement fields is here stud- 
ied in the harmonic approximation by introducing the 
Fourier series u(r) = 2 Uqe*'! '" and h{r) — 

'^qS*''''^! where V is the volume of the system. 



A. In-plane phonons 



The decoupled in-plane phonon modes are obtained in 
the usual way by changing to longitudinal Uq = Uq • q/g 
and transverse Uq = Uq • (ez x q/f?) displacement fields. 
The dispersion relations have the usual linear behavior 
in momentum and are given by 



Wq = VLQ, 
Wq = VtQ, 



(3) 



2. Strained case 

Suspended samples may be under tension either due 
to the load imposed by the back gate or as a result of 
the fabrication process, or both. The case of a clamped 
graphene membrane hanging over a tranche of size L, 
relevant for conventional two-contacts measurements in 
suspended samples, has been considered in Ref. [l^. 

Once the membrane is under tension a static defor- 
mation configuration is expected at equilibrium. The 
phonon modes may be obtained by assuming that both 
in-plane u(r) and flexural h{r) fields have dynamic com- 
ponents which add to their static background: u(r) ~ 
(r) and h{r) = hst{r) + hdyn{r). For the 
case of the clamped membrane considered in Ref. [l^ we 
have Ust(r) = [ux^stix),0] with Ux^st(x) a linear function 
of X, while hst{x) may be approximated by a parabola. 

Let us consider the general static displacement vector 
field dst(r) = [u^^st(r), ^^^^^((r), /ist(r)] and the associated 
strain tensor u^j,st{v) = \{diUj^st + djU^^st + dihstdjhst)- 
In-plane phonons are not affected by the static compo- 
nent but the FP dispersion changes considerably. This is 
a consequence of new harmonic terms appearing due to 
coupling between in-plane static deformation and out-of- 
plane vibrations in the full strain tensor in Eq. ([2]). The 
resulting FP dispersion relation may be obtained using a 
local approximation expected to hold for L ^ £ = vft ^ 
kp^, where vp is the Fermi velocity, kp the Fermi mo- 
mentum, and T a characteristic collision time. The result 
reads 



A 



\ p p p 



(5) 



For the particular case of isotropic strain where Ux 



'■yy 



and Uxy = the dispersion relation can be cast in 



the form 



t^q(i") = q\ -q'^ + uu,st{r) 



(6) 



with vl 



and vt 



Typical values for 



monolayer and bilayer graphene are vl — 2.1 x 10*m/s 
and VT — 1.4 x lO^m/s. 



B. Flexural phonons 

1. Non-strained case 

The quadratic behavior in out-of-plane displacements 
of the strain tensor in Eq. ([2]) implies that FP modes are 
driven by the bending rigidity term. The resulting FP 
dispersion relation is quadratic jii^i^ 



aq 



(4) 



Here we will give particular emphasis to the clamped 
membrane case where the FP dispersion is given by 



A -I- 2^ _2fi , 2 
u — sm I 



P 



P 



(7) 



with u = Uxx and (/)q ~ a,TCia,n{qy / qx) ■ In order to keep 
the problem within analytical treatment we will use an 
effective isotropic dispersion relation, obtained by drop- 
ping the angular dependence contribution. 



qya-q^ + uvl- 



(8) 



Since we are mainly interested in transport such an ap- 
proximation has the advantage that backward scattering 
is still correctly accounted for. 
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III. ELECTRON-PHONON INTERACTION 



A. Low energy description for bilayer graphene 



indices and reads 

^i = E^k [^i(k,k')I]Vk', 



(12) 



At low energies the 2-band effective model provides a 
good approximate description for tt— electrons in bilayer 
graphenciii The 2x2 Hamiltonian can be cast in the 
form Ueff = I]k V'k'HkV'k- with 







2ni \ (kj. + ikyY 



(kx iky)' 




(9) 



where the two component spinor ipj^ = [aj^, stems from 
the two sublattices not connected by the interlayer hop- 
ping, t± « 0.3 eV. The coupling t± between layers sets 
the effective mass 2m = t±/vp, with vp ~ lO^m/s the 
Fermi velocity in monolayer graphene. Equation ^ is 
valid at valley K, at valley K' we have — >■ 'Hj ■ Here 
we are interested in electron scattering processes induced 
by emission or absorption of long wavelength acoustic 
phonons and hence intervalley scattering is not allowed. 
Thus we may concentrate on one valley only. The Hamil- 
tonian can be diagonalized introducing the rotated oper- 
ators 



(10) 



1 / fo'Sk p~i^k 

where 0k = arctan(fcj,/fca;), with 
such that ejj. stands for electron-like (positive energy) ex- 
citations and /ij^ for hole-like (negative energy) excita- 
tions, from which we get 



[ejj., hl^] defined 



n 



k 



e (k) [e^ ek ■ 



(11) 



with e(k) = h^k'^/{2ra). 

Along the paper we will compare the results ob- 
tained for bilayer graphene with those valid in mono- 
layer. The latter is described using the effective Dirac- 
like Hamiltonian^ii which holds in the low energy sector 
we are interested here. 



B. Coupling between electrons and phonons 

The coupling between electrons and the vibrations 
of the underlying lattice either in bilayer or monolayer 
graphene has two main sources. Long wavelength acous- 
tic phonons induce an effective local potential called de- 
formation potential and proportional to the local con- 
traction or dilation of the latticcjii 

^i(r) = gouuiy), 

where is the bare deformation potential constant, 
whose value is in the range go ~ 20 — 30 eV.^" The re- 
spective interaction Hamiltonian is diagonal in sublattice 



where I is the 2x2 identity matrix and Vi(k, k') is the 
Fourier transform of the deformation potential 



Fi(k,k')=V-\9o j dve'^'^'-'^y^uu (r) 



(13) 



Equation (jl2p is valid both for monolayer and bilayer 
graphene. Since we are interested in doped systems we 
take into account screening by substituting Vi(k, k') with 
Vi (k, k') /e (k — k'), where we take a Thomas-Fermi like 
dielectric function 



e(q) = i 



2eoq 



(14) 



and V [Ep) is the density of states at the Fermi energy, 
which is given by V (Ef) = -M^ = in the case of 

in the case 



monolayer graphene and by V (Ep) 
of bilayer. 

It is convenient to define g = go/e{kF) for single layer 
graphene, which gives a density independent screened de- 
formation potential 



9 ■■ 



9o 



9o 



e^/ineohvp) 8.75 



2 - 3.5 eV. 



(15) 



Note that the value just obtained is in complete agree- 
ment with recent ab initio calculations which give g « 
3 eVj^ It will become clear in Sec. IIV C II that g as de- 
fined in Eq. ([T5|) is the relevant deformation potential 
electron-phonon parameter in single layer graphene. For 
bilayer graphene g(q) = go/e{q) gives 



9ikF) 



90 



90 



eH±/{2Treoh?vj,kp) 11.25 



a (2-3)^^: eV, 

with n in 10^^ cm ^. We may then write a q dependent 
deformation potential electron-phonon parameter which 
has the form gMil) = 9<l/kp for monolayer graphene, 
and gsil) = 92hvpq/t± for bilayer. 

Phonons can also couple to electrons in monolayer and 
bilayer graphene by changes in bond length and bond 
angle between carbon atoms. In this case the electron- 
phonon interaction can be written as due to an effective 
gauge field r^'"— 



eA, 



elastic 



hiux 



where /3 « —dlogt/dloga ^ 2 — 3^ with t the in-plane 
nearest neighbor hopping parameter and a the carbon- 
carbon distance {to sa 3 eV and oq ~ 1.4 A). In the case of 
bilayer graphene the resulting interaction Hamiltonian is 
obtained by introducing the gauge potential into Eq Q , 
following the minimal coupling prescription, and keeping 
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only first order terms in electron-phonon coupling. Then 
we arrive at 



2m ^ ^ 

k.k' 







k.k' 



k.k' 



(17) 

where tt^ = A:e±'''^ A^^, = y2,.(k, k') ±iF2,y(k, k') and 
the vector V2 = (V2,x, 1^2, y) is defined as 

y2,.(k,k') =^V-i y dre^e^'-'')-^^ [«,,(r) - u^^(r)] , 

y2,,(k, k') = ~ ^V-i y dre'('''-'')-u,,(r). (18) 

An estimate of the electron-phonon coupling strength due 
to V2 is given by kpfi^fi/irna) sa (8 — 12)Y^eV, with n 
in lO^^cm"^. In the case of single layer graphene the 
resulting interaction Hamiltonian reads 



H2 =l^F^V^lT.V2(k,k')^k' 



(19) 



k.k' 



where cr = {<7x,ay) is the vector of Pauli matrices, and 
the two component spinor -0^ = [aj^., is reminiscent 
of the two sublattices of the honeycomb lattice. An esti- 
mate of the respective electron-phonon coupling strength 
is given by vpfiP/a w 10 — 15 eV. 



The electron-phonon interaction Hamiltonian is the 
sum of the two terms shown above, iJep = Hi -|- ^2- 
Phonons enter through the strain tensor Uij which we 
have seen can be written in terms of static and dynamic 
components; the static ones being zero for zero load. 
There are purely static terms which do not contribute 
to electron-phonon scattering and will be dropped (see 
Ref. [Tgl ). Quantizing the dynamic part of the displace- 
ment fields;^ and introducing usual destruction and cre- 
ation operators and (flq)^ for in-plane phonons q and 
polarization v = L,T, we can write the q component of 
the in-plane displacement as, 



L/T 



L/T 



(20) 



For FPs we introduce the bosonic fields and [a^)\ 
and write the q component of the out-of-plane displace- 
ment as, 



/iq = 



(21) 



The electron-phonon interaction Hamiltonian may then 
be written either in monolayer or bilayer graphene as. 



Hp' 



k,k' [ f,q q,q' J 



E \ E "^^Ab^' K + («-q)1 '^k'.k-q + E "^2,^,^^^' [< + («'^q)1 [<' + (^'^qOl -^k'.k-q-q' + h.C. 

k,k' I y,<i q,q' 



(22) 



For monolayer graphene the matrix elements read, with = (see also Refs. Pzlfisl). and where we have 

again used the local approximation. In the case of bi- 
layer graphene only the matrix elements for the gauge 



^^'^ ~ efo) *'^\/ 2V ouj^ ' potential change, becoming dependent on fermionic mo- 

menta k and k'. As can be seen by comparing Eqs. P7)) 
\ ^ and ([T^. they take exactly the same form as in Eq. 




AV p^l Lo^ uj^i with the replacement vp — > fi(7rk + Tiy^i) / {2m) . 



mF 90 . ^ , 



1 e{q) y 2Vpuj^' IV. TEMPERATURE DEPENDENT 
RESISTIVITY 



^'2,q - ^„ 



2a V/ 2Vpa;^ Our aim here is to study the T dependent resistivity 

I in suspended bilayer graphene as a result of the electron- 

yT _ ^^-^^ g(^'20a / ^ ^ phonon interaction derived above. We assume the doped 

2a Y 2VpWq regime iSp ^ /i/r, where l/r is the characteristic elec- 

p _ , ,;(0 



4a ov^ /,.,F,.,F 



2Vp./<^^ 
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tronic scattering rate (due to phonons, disorder, etc). 
The doped regime iniinediately iniplies kp^ <^ vpT = 
where / is the characteristic mean-free path, thus justify- 
ing the use of Boltzmann transport theory (even though 
graphene's quasiparticles are chiral the semiclassical ap- 
proach still holds away from the Dirac pointi22i2^) . 



A. The variational approach 

The Boltzmann equation is an integro-deferential equa- 
tion for the steady state probability distribution /k>^ It 
can be generally written as 



and 

($,P$)= *k7'k,ki,...,k„($k±3'ki---±*kj 

k.ki,...,k„ 



^ ^ ($k±1>k,---±$kj'Pk,k,....,k„ 

(29) 



(n+1) 



k,ki,...,k„ 



the variational principle asserts that of all functions 
(f>k satisfying = ((f>,P$), the solution of the 

linearized Boltzmann equation gives to the quantity 
its minimum valuej^ In particular, 
the resistivity g can be written as 



i- • Vr/k + k • Vk/k = /k 



(24) 



1 



where the terms on the left hand side are due to, respec- 
tively, diffusion and external fields^ while on the right 
hand side scattering provides the required balance at the 
steady state. (See Appendix |X] for an explicit form of 



/k 



in the case under study.) The Boltzmann equa- 



tion is quite intractable in practice, and its linearized 
version is used instead, 



k • Vk/, 



(0) 



5h 



(25) 



where (^/klg^^^^j. is the linearized collision integral. (See 
Appendix |B] for an explicit form of 5f\A in the 

' I SC3.tt 

case under study.) Expanding the distribution proba- 
bility around its equilibrium value f^'^ ~ l/{cxp[(ek — 
fl)/kBT] + 1}, 



AO) 

Jk 



(26) 



0. it 



and using the equilibrium property that /A i 
can be seen that S f-kl is linear in <f>k. and that it can 

^ Iscatt ... 

be written as a linear application in terms of the linear 
scattering operator Pk, 

5/k|_,, = Pk1>k^ 

= - J2 ^k,ki,...,k„ (*k±*ki •••±*kj, 

ki....,k„ 

(27) 

where 7^k,ki,...,k„ is a generalized transition rate per 
unit energyi^ (See Appendix [B] for an explicit form of 
'Pk,ki,...,k„ in the case under study.) Writing the lin- 
earized Boltzmann equation, Eq. (|25p . in the form 



Xk 



k^k, 



and defining the inner products, 
(<i>,X)-^$k (r 



k • Vk/r 



(28) 



gd{($,X(i? = l,Vr/(o) = 0))}2' 



(30) 



being thus expected to be a minimum for the right 
solution;^ where gd is the system's degeneracy (for 
monolayer and bilayer graphene it is g^ = g^g^ — 4 
due to spin and valley degeneracies). The quantity 
X{E = l,Vr/(°) = 0) refers to the left hand side of 
Eq. ((25| in a unit electric field and no spatial gradients 
(for example, zero temperature gradient). It is easy to 
show that ($, X) = E • J, where 



J = ^ evk$k 



dew 



is the current per non-degenerate mode (per spin and 
valley in monolayer and bilayer graphene). The quan- 
tity {{<i>,X{E = 1, Vr/(°) = 0))}2 is therefore nothing 
but Vj^, where j — J /V is the current density. 

A well known solution to the Boltzmann equation ex- 
ists when scattering is elastic, the Fermi surface isotropic, 
and the transition rate can be written as Pk,k' = 
'P(fc, ^^k.k')i where 6'k.k' = ^k — ^k' is the angle between k 
and k'l^ Under these conditions the solution reads, 

<i>k = I'k • (eE - ^vr) T{k), 

where T{k) is the isotropic scattering rate, and we have 
written Vr/^"^ = df'^^/dei^VT. Clearly, the later solu- 
tion for $k can be cast in the form <I>k (x k- Uf^i^S where 
u is a unit vector in the direction of the applied fields. 
So, in more complicated cases where there is a departure 
from the isotropic conditions and/or from elastic scat- 
tering, it is a good starting point to use Eq. ([50]) with 
$k c)c k • u to get an approximate (from above) result 
for the resistivity. Note that the coefficient multiplying 
k • u is unimportant as it cancels out. This variational 
method is equivalent to a linear response Kubo-Nakano- 
Mori approach with the perturbation inducing scattering 
treated in the Born approximation.— 

Here we use the variational method just outlined to 
get the T dependent resistivity in bilayer graphene (and 
monolayer for comparison) due to scattering by acoustic 



6 



phonons. In this case, using the quasi-elastic approxima- 



tion (see Appendix [B]), (5/k| 
form of Eq. ^ 



can indeed be cast in the 



5h\ 



(31) 



where for scattering by one in-plane phonon 



Pk,k' = 



Wi.(q, k,k )t^q^-r^ ()k,k'+qO(ek-ek' 



and for scattering by two FPs 
2n. „a/f 



(32) 



T'k 



knT- 



dei. 



(q, q , k, k ) ^^-^ -:— X 



F F 
^q - ^q' 

nq - 71q' 



^k,k'+q+q"5(ek — £k' 



(33) 



with nq = l/[cxp(?ia;q/fcBT) — 1] the equilibrium phonon 
distribution. The kernel quantities it;i/(q, k, k') and 
^^-(q, q', k, k') are related to the matrix elements in 
Eq. as follows (see Appendix for bilayer 

graphene, 

«;,(q,k,k') « iFi-^ql' (1 + cos20k,k') 

2 



2,q 



fc'' 



2fcfc'cos6'k,k') . (34) 



for one phonon processes, and a similar expression for 
two phonon processes Wi?(q, q', k, k') with — )■ V^^,, 

where V2 means the matrix elements given in Eq. (|23p 
for the gauge potential without the term (vTk -I- TTk'); for 
monolayer graphene, in the case of one phonon process. 



u;,(q,k,k') « iFi-^ql' (1 + cos^k.k') + IV^a^ql' 



(35) 



with a similar expression for two phonon processes 
«;f(q,q',k,k') withVq^'^Vq^q,. 

Using the setting given above the resistivity is conve- 
niently written as 

1 5Ek,k'(*k-$k')'a,k' 



5s 



V 



V /dkdk'(K-u)2-Pk,k' 



J dk k • uvk 



(36) 



where we changed from summation over k— space to 
integration, and defined K = k — k'. The integral 
in the denominator can be done immediately assuming 
ep 2> ksT. The result reads the same for bilayer and 
monolayer graphene. 



/ 



dk k • uVk 



(37) 



In order to proceed analytically with the integral in the 
numerator we have to specify the T regime, as discussed 
in the next section. 



B. Bloch -Griineisen temperature 

For each scattering process (one or two phonon scat- 
tering) we may identify two different T regimes, low and 
high T, related to whether only small angle or every an- 
gle are available to scatter from |k) to |k'). Recall that 
since we are dealing with quasielastic scattering both k 
and k' sit on the Fermi circle, see Fig.IB and |k) and |k') 
are adiabatically connected through a rotation of ^^k.k' in 
momentum space. Large angle scattering is only possi- 
ble if phonons with high enough momentum are available 
to scatter electrons. The characteristic Bloch-Griineisen 
temperature Tbg separating the two regimes may thus 
be set by the minimum phonon energy necessary to have 
full back scattering. 



'2kf 



(38) 



with ujq as given in Sec. HIl 

For scattering by in-plane phonons Tbg takes the value 



T 



(L) 
BG 



57V^K and t'^' « 38V^K 



(39) 



for longitudinal and transverse phonons respectively, 
with density n in units of lO^^cm"^. When scattering 
is by two non-strained FPs, the crossover between low 
and high T regimes is given by 



Tbg «0.1nK, 



(40) 



with n again measured in lO^^cm"^, while in the pres- 
ence of strain, using the approximated strained FP dis- 
persion in Eq. ([8]), we get 



Tbg — 28vmi K. 



(41) 



It is obvious from Eqs. (gOl) and (gT]) that the high-T 
regime is the relevant one for FP scattering. 



C. Contributions to resistivity 

In the following we summarize our results for the T de- 
pendent resistivity due to scattering by in-plane phonons 
and two FPs in bilayer graphene. For comparison we dis- 
cuss also the monolayer case first studied in Ref. 0- We 
use the variational method discussed in Sec. lIV Al the re- 
sistivity being given by Eq. p6p . Details on the deriva- 
tion can be found in Appendix [Cj We neglect one FP 
processes since these, as can be seen in Eq. (|23p . are re- 
duced by a factor ^ Hq/L <^ 1, where Hq is the sample's 
vertical deflection over the typical linear size L. 
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q' 



where T{n) = (n — 1)! is the gamma function and (^{n) 
is the Riemann zeta function. We have thus obtained 
the expected behavior at low T for couphng through 
gauge potential, which is the 2-diniensional analogue of 
the Bloch theory in S-dimensional metals»^i^ The 
scalar potential contribution comes proportional to 
due to screening. It can be neglected in the low T regime; 
even though 16r(6)C(6)/[r(4)C(4)] w 300, it is strongly 
suppressed by T/Tbg ^ 1 and g < hvFl3/{2a), (see 
Sec.MEl- 

In the high T regime, T 3> Tbq, the inequality <C 1 
holds, so that 6^-^/(6^-^-1)2 ^ l/{z^x) inEq. The 
usual linear in T resistivity for one phonon scattering is 
then recovered, 



Figure 1: Scattering of electrons in momentum space due 
to: (a) in-plane phonons, (b) non-strained flexural and 
(c) strained flexural phonons. 



1. Scattering by in-plane pfionons 

A sketch of the scattering process in momentum space 
involving one phonon is shown in Fig. [IJa) . In this case 
the resistivity can be written as (see Appendix I C ip . 



8hkl 



1 



(42) 

where ^ huj2i^^/kBT, with v = L,T, and where we 
have introduced a generalized electron-in-plane phonon 
coupling for bilayer graphene given by. 



29V 



4 



-,1/2 



(43) 

The case of scattering via screened scalar potential, here 
encoded in the screened deformation potential parame- 
ter g, has been considered recently in Ref. |32]. As it is 
shown below, the gauge potential contribution becomes 
the dominant one in the low T regime. 

In the low T regime, T ^ Tbg, we have 3> 1, so that 
the integrand in Eq. (|42p is only contributing significantly 
for a; ^ 1. The generalized electron-in-plane phonon 
coupling in (|33]) then becomes. 



^b(2/«1) 



1/2 



(44) 



and the resistivity reads. 



E 



,16r(6)C(6) / T 



r(4)C(4) \T, 



'-BG 



4a2 



r(4)c(4)(fcBr)^ 



(45) 



nkBT 



(46) 



where l/w^ = + Note that, at odds with 

the low T regime, now the scalar potential contribution 
is higher than the gauge potential one for the typical 
coupling values discussed in Sec. IIII Bl 

The monolayer case has been discussed extensively in 
the literature.— liii^i^i^i The resistivity is still given by 
Eq. (|42p . only the generalized electron-in-plane phonon 
coupling changes. 



29V 



1 



r 
4 



4a2 



1/2 



(47) 



The same qualitative behavior is obtained: at low T the 
resistivity is given by Eq. (|45p with the numerical re- 
placement 16 — >■ 12 in the scalar potential contribution; 
at high T the result (j^Hl) holds with the replacement 

+ - + Note, how- 

ever, the apparent quantitative difference: the scalar and 
gauge potential contributions change roles, the later be- 
coming more important in monolayer graphene. This is 
further discussed in Sec. |Vl 

A final remark regarding the temperature dependent 
resistivity due to in-plane phonons has to do with the 
value of the electron-phonon coupling parameters /3 and 
g. While (3 is expected to be restricted to the range 
P ^ 2 — Z, as discussed in Sec. IIII Bl the value of the 
deformation potential parameter g is still debated in the 
literature. Phenomenology gives g ^ 10 — 30 eVj^i^ re- 
cent ab initio calculations provide a much smaller value 
gr ~ 3 eVi^ On the other hand, experiments seem to 
confirm the higher values, giving g '--^ 15 — 25 eSf 
Our claim here is that all these values make sense, if 
properly interpreted: phenomenology gives essentially 
unscreened deformation potential, which we called go in 
Sec. IIIIBl and which should take values of 0(10) eV; 
screening effects suppress the deformation potential to 
0(1) eV, as we have seen in Sec. HiTB] within the Thomas 
Fermi approximation, in good agreement with ah initio 
results where screening is built in; the fact that trans- 
port experiments give a much higher deformation poten- 
tial is a strong indication that phonon scattering through 
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Figure 2: Resistivity vs T due to scattering by in-plane 
phonons (dashed thick lines) and flexural phonons (full lines) 
with and without strain (indicated in percentage), in mono- 
layer (a) and bilayer graphene (b). 
g K. 3eV, and /3 3. 



We use n = 10 cm 



gauge potential, usually not included when fitting the 
data^i^^ is at work. Indeed, using the monolayer ver- 
sion of Eg. P5|) , we readily find that the fitting quantity 
in Refs. |35| and|3^ should be replaced by. 



D 



2a2 



1/2 



pling. 



Dm(x) 



g'x' ( 1 



X 

T 



4a2 



1/2 



(51) 



In Fig. [T]^b) a sketch of the two phonon scattering pro- 
cess in momentum space is provided. It shows that one 
of the two phonons involved in the scattering event al- 
ways has momentum q' — > 0. This is a consequence 
of the quadratic FP dispersion [Eq. ([4])], which leads to 
a divergent number of FPs with momentum q' — >■ 0.— 
This divergence is responsible for the logarithmic factor 
in Eq. (j49p , which stems from the existence of an infrared 
cutoff qc- This cutoff is to be identified with the onset of 
anharmonic effects,— or unavoidable built in strain;^ 

In the low T regime, T ^ Tbg, one has z ^ 1, so that 
the integrand in Eq. (j49|) is only contributing for a; ^ 1. 
The generalized electron-phonon coupling becomes equal 
in both bilayer and monolayer systems. 



■ 2 2 , fl'vll3^''^' 

9 y ' 



4a2 



(48) and the resistivity is then the same in both. 



(52) 



which, keeping g ~ 3 eV, takes values D ~ 10 — 20 eV 
for /? ~ 2 — 3, in excellent agreement with experiments. 
Moreover, since the gauge potential is not screened 
[Eq. it provides a natural explanation for the 

resistivity beliavior recently reported at low T in Ref. IsgI . 
where the expected contribution due to scalar poten- 
tial is absent 1^ 

Typical T dependent resistivity due to scattering by 
in-plane phonons, Eq. (|42|) . is shown in Fig. [2] as thick 
dashed lines. In agreement with analytical results, there 
is no qualitative difference between monolayer graphene 
[Fig.ll^a)] and bilayer [Fig. [D^b)]. 



9F 



■ 6r(6)c(6) 
r(4)c(4) 

r(4)c(4);ifc 



T 



T, 



BG 



hakp 



4a2 

5/2 

In 



ksT 
haq^ 



(53) 



A similar result has been derived in Ref. 13 Owing to 
the same arguments used in the previous section for one 
phonon scattering we can neglect the scalar potential con- 
tribution at low T . 

At high T, i.e. T > Tbg, we have z < 1, so that 
exp(za;2)/[exp(za;2) - 1]^ « l/{zx'^) in Eq. The 
bilayer graphene resistivity becomes. 



2. Scattering by non-strained flexural phonons 

In non-strained bilayer graphene scattering by FPs give 
rise to the following T dependent resistivity (details on 
the derivation are given in Appendix I C 2|) . 



In 



ksT 
haq^ 



dx ^ 
Vl 



a;2 (e^ 



where z ~ fujj2kp/kBT, and where the generalized 
electron-FP coupling for bilayer graphene is given by, 

-I 1/2 



D'b{x) 



2 2 

9 X 



1_ ' 



4a2 



1-^ 

4 



(50) 

Equation (|49p holds also for monolayer grapliene, we need 
only to introduce a different generalized electron-FP cou- 



QF 



{kBTf 



QAhe^p^vla^kj, 



In 



kBT 
haq^ 



(54) 

This result holds for monolayer graphene with the sub- 



stitution g 



^ -1)2' 

(49) 



Ba^ J ' ' ^j^'iWehave 
obtained that the resistivity due to non-strained FPs 
is proportional to T"^ /n, which implies mobility inde- 
pendent of the carrier density n. A similar result has 
been obtained in the context of microscopic ripples in 
graphenci^i^ The result of Eq. ((M)) is shown in Fig. [5Jb) 
as a full line indicating u k, 0%. The logarithmic cor- 
rection, expected to be of order unity in the relevant T 
range, has been ignored. Scattering by FPs dominates 
the contribution to resistivity in non-strained samples 
at both low and high T, except for the crossover region 
where T ^ Tbgi Eq. ((5^ . The same conclusion holds for 
monolayer graphene, whose T dependent gp is shown in 
Fig.UTa). 
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3. Scattering by strained flexural phonons 

Applying strain breaks the membrane rotational sym- 
metry inducing linear FP dispersion at low momentum, 
as can be seen in Eq. ([S]) . A new energy scale appears in 
the problem, 

wf. = V2uvl/a « 10*u(K), (55) 

separating two regimes: linear dispersion below and 
quadratic dispersion above. The associated momentum 
scale, q* = y/UvL/a ~ 4.5-\/s A , together with kp and 
the thermal momentum given by = k^T, define 

all regimes where analytic treatment can be employed. 
In particular, in the low T regime where qx ^ kp we 



may always take q* ^ qx and use a linear dispersion 
for FPs; otherwise the non strained case considered in 
the previous section would be the appropriate starting 
point. In the high T regime we can distinguish between 
low strain for q* <C qr and high strain for q* ^ qr- Note 
that at high T relevant phonons scattering electrons have 
momentum q in the range kp < q < qx- Therefore, when 
strain is present in the high T regime we may always 
assume q* ^ kp] the opposite limit, q* <^ kp, would 
again be identified with the non-strained case considered 
previously. 

The resistivity due to strained FPs can be cast in the 
form of a triple integral over rescaled momenta a; — )■ i = 
xfiVLU^/"^ /{kBT) (see Appendix IC 21 for details). 



Qf 



{kBTf 



2^'K'^h^e'^p'^v1,v\u'^k\ jQ 



dK'- 



dq —nq{nq + l)x 



\K+q\ 



dQ- 



K-q\ 



ujA\lfK^- [K^ + q^-Q^] /4 

r 



1 + ng + Uq 



Tig - Uq 



(56) 



where the rescaled dispersion reads utx « j^x^ + i^, 
with 7 = \/2L0g^ I ujg, ^ and the generalized electron-FP 
coupling Dg{y) is given by Eq. (j50p . For monolayer 
graphene only the coupling changes, being given instead 
by Eq. ((5T|l . The kinematics of the scattering process is 
schematically shown in Fig. [IJc) . 

In the low T case, T <C Tbg, we have only small angle 
scattering with K <^ kp. The argument of the gener- 
alized electron-FP coupling becomes small, K/kp ^ 1, 
and it can be written as 



^s(?/< 1) «i^A/(y< 1) 



2 2 

9 V 



4a2 



(57) 

The resistivity is the same in both bilayer and monolayer 
systems. Since the inequality qp <^ kp,q* holds, rele- 
vant phonons have linear dispersion ojg w \/Uvpq and 

the rescaled Fermi momentum obeys kp kp/qp 3> 1. 
We may take K ^ oo as the upper limit in the K integral 
in Eq. (|56p . and the resistivity is then approximated by. 



where 



ICr, 



dKK" 



dqq nq{nq + l)x 



dQ- 

k-q\ 



q+Q 



q-Q 



q^K^ — 



2 q2 



/4 



(59) 



It can be shown numerically that IC2 ~ 4485 and /C4 w 
496850. The large ratio K.4/K.2 ^ 1 is, however, com- 
pensated hy qp/kp <^ 1 and the fact that g < hvpf3/{2a) 
(see Sec. IIII B|) . As in the case of scattering by in-plane 
phonons, also here the gauge potential contribution to 
resistivity dominates at low T. 

Now we consider the high T regime, T ^ Tbg- At 
odds with the non-strained case [see Fig. mb)], now 
phonons with momentum q in the range kp < q < qp pro- 



vide most of the scattering. It is shown in Appendix IC 2 bl 
that the integral over Q in Eq. ([55)1 becomes K indepen- 
dent, and the q integral can be cast in the form. 



QF 



51 



4a2 



{kBT)"' 



dq 



^ q'niinq + 1)^ 



g2 



1 



1 + 2nr, 



nq{nq + 1) 

267r2n6g2p2^|y9 ^9/2^3^' ~- ^gp-j 

(58) being easily evaluated numerically. The resistivity in bi- 
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layer graphene can then be written as 



QF 



7.r 



4a2 



2'^ p^v^v^i 



huv\ 



(61) 

When 7 ^ 1, or equivalently qT ^ q* (high strain), the 
function ^(7) behaves as ^(7 < 1) « 18C(3) - 93C(5)/8. 
For 7 ^ 1, or equivalently qx 3> q* (small strain), it 
gives t/(7 S> 1) ~ 1/7^. In these asymptotic regimes 
one can obtain analytic expressions for the resistivity in 
Eq. 161]), 



QF 



1 



4a2 / 27ft3e2„2 



f(l8C(3)-fC(5))5|fg kF^qr^q* 



kp <^ q"" <^ qt 



(62) 



Equations (pT|) and (|62|) also hold for monolayer graphene 



with 



2.9^ 



The effect of strain in the T dependence of resis- 
tivity is shown in Fig. [H^a) for monolayer graphene 
and[5][b) for bilayer at strain values u « 0.1%, 0.5%, 1%. 
The crossover between the two regimes of Eq. (pT]) [see 
Eq. (HH)] is clearly seen at 7 sa 1, which is equivalent to 
T m lO^u K. It is apparent from Fig. [2] that the con- 
tribution to the resistivity due to scattering by FPs is 
strongly suppressed by applying strain. 



D. Crossover between in-plane and flexural phonon 
dominated scattering 

Scattering by in-plane and flexural phonons are always 
at work simultaneously. However, the two mechanisms 
provide completely different T dependent resistivity, and 
therefore we expect them to dominate at different T. In 
the following we address the transition T at which ~ 



1. Non-strained case 

In this case, using Eqs. and ([M]) for bilayer 

graphene in the high T regime we get 



QF 



T[K] 



Qm 75n[cm 2] ' 



(63) 



We expect a crossover between in-plane to FP dominated 
scattering given by 



Tc2 « 7571 [cm"2] X. 



(64) 



The Tc just obtained is close to Tbg for in-plane phonons, 
Eq. dMl)- and much higher than Tbg for FPs, Eq. pi)l . 



Using the low T approximation for Eq. (|45p . we ob- 
tain the ratio 



QF _ , , „ y^^pri^ 



12- 



Qv 



{T[K]f 



(65) 



from which we expect a crossover from FP to in-plane 
dominated scattering at 



T,i « 3(n[cm-2])i/4 k. 



(66) 



as T increases. We conclude that scattering by FP al- 
ways dominates over scattering by in-plane ones, except 
for the region Td <C T <C Tc2 around Tbg for in-plane 
phonons, Eq. ([M]) . This is clearly seen in Fig. HKb). 
The same conclusion applies to monolayer graphene. In 
this later case we obtain Td ~ 6(n[cm~2])i/4 k and 
Tc2 « 55n[cm-2] K. 



2. Strained case 

It can easily be shown that the crossover from in-plane 
to flexural phonon dominated scattering always occurs 
in the low strain regime, q* <^ qT- We have seen in 
the previous sections that the crossover temperature T 
separating high strain from low strain behavior is given 
by 7 = \f2ijj^^ /u^, « 1. Using Eq. we get a crossover 
temperature T* k, 10'*uK. On the other hand, using 
the low strain approximation for the resistivity due to 
flexural phonons given in Eq. (j62p and the resistivity due 
to in-plane ones in Eq. (|46p we obtain for the ratio in 
bilayer graphene 



QF _ ksT 

Qin AOtTKU 



The corresponding crossover T then reads 



Tr w 10^1 K. 



(67) 



(68) 



Clearly Tc 3> T*, justifying our low strain approxima- 
tion. The same applies to monolayer graphene under 
strain. The resistivity ratio is in that case QF/Qin ~ 
kBT / [bQ'KK.u), from which we obtain roughly the same 
Tc even taken into account that k in monolayer is half 
that of bilayer in our elasticity model. 

In important conclusion may be drawn. While in the 
non-strained case scattering by FP is the dominant con- 
tribution to the resistivity, it can be seen from Eq. (j68p 
that applying small amounts of strain is enough to sup- 
press this contribution at room T . This is clearly seen in 
Fig.H 



V. DISCUSSION 

We have found the T dependent resistivity due to 
acoustic phonons to be qualitatively similar in monolayer 
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and bilayer graphene (see Sec. IIV C[) . This becomes ap- 
parent when we compare Fig.[5Ja) andl^fb) where g{T) is 
shown, respectively, for monolayer and bilayer graphene 
both at zero and finite strain. Such behavior can be 
traced back to the resistivity expression, Eq. ([55)1 . or 
more precisely to its numerator, where the different elec- 
tronic structure and electron-phonon coupling conspire to 
give exactly the same parametric dependence in mono- 
layer and bilayer graphene. In short, it can be read- 
ily seen through Eq. ([M]) that the information about 
the electronic structure enters via the density of states 
squared in the numerator. The electron-phonon cou- 
pling, in its turn, enters through the transition rate Vk^k' , 
and in the Born approximation it also appears squared. 
When coupling is via scalar potential [Vi matrix elements 
in Eq. (|23p] the screening makes the coupling inversely 
proportional to the density of states, and the square of 
it cancels exactly with the square density of states com- 
ing from the integral over k and k' in the numerator of 
Eq. ([55)1 . For the coupling through gauge potential [V2 
matrix elements in Eq. ()23p ] the parametric difference 
between monolayer and bilayer amounts to the replace- 
ment Vp — >■ h'^kp/{2m)'^ (after taking the square of the 
matrix elements and using the quasielastic approxima- 
tion). When multiplied by the square of the density of 
states in bilayer graphene, I?(i?) ^ m/fi^, we obtain the 
factor kp/h?, which has exactly the same parametric de- 
pendence appearing for single layer graphene - there, the 
factor Vp multiplies the square of the density of states of 
the monolayer, !?(£') ~ k/(hvp). 

Despite qualitative similarities there are apparent 
quantitative differences. A striking one is the overall 
suppression of resistivity in bilayer graphene, which is 
clearly seen when we compare Fig. |2][a) and lljb) . This 
is due to the higher stiffness and mass density of bilayer 
graphene (see Sec. |ll|), and to the 1/2 term in Eq. ([TT)) 
for the gauge potential which, at odds with the paramet- 
ric dependence just discussed, does not cancel out in the 
expression for the resistivity. Also, the scalar potential 
contribution is quantitative different, being enhanced in 
bilayer graphene, as can be seen by comparing Eqs. ()45|l . 
and dM)), Eqs. ^ and dMj), and Eqs. §^ and §^ 
with their monolayer counterparts. This arises because 
pseudo-spin conservation allows back scattering due to 
scalar potential in bilayer but not in monolayer graphene. 

The quantitative discrepancy between the T depen- 
dent resistivity in bilayer and monolayer graphene origi- 
nates an interesting difference regarding room T mobil- 
ity in non-strained samples. The mobility /i, defined as 
g = l/(en^), is in the non-strained case limited by PF 
scattering (see Sec. IIV C 2)1 and takes the form 



GAnhevp 



(69) 



with Ab = IW + l{%a'')\ in bilayer graphene 

and Am = '«^/[3^/2 + fi^Vpfi"^ /{Aa"^)] in monolayer 
graphene, ignoring the logarithmic contribution of or- 
der unity, Eq. ()54p . For monolayer graphene at room 



T the mobility is limited to the value for samples on sub- 
strate, fi « 198Am ~ 1 m^/Vs, as has recently been con- 
firmed experimentally.— For bilayer graphene, however, 
the quantitative differences discussed above lead to an en- 
hanced room T mobihty, ^ « 198Ab ~ 20 m^/Vs. This 
might be an interesting aspect to take into account re- 
garding room T electronic applications. Reports of much 
smaller mobility (one order of magnitude) in recent ex- 
periments in suspended bilayer graphenei might be an 
indication that residual, T independent scattering is at 
work, overcoming the intrinsic FP contribution. 

Another remark worth discussion is the validity of re- 
sults in the non-strained regime, in particular Eq. (|54p 
for the high T resistivity. At low densities, when kp be- 
comes comparable with the infrared cutoff Qc given by 
the onset of anharnionic effects,— the harmonic approxi- 
mation used here breaks down. A complete theory would 
require taking into account anharmonicities, but this is 
beyond the scope of the present work. Nevertheless, it is 
likely that unavoidable little strain u is always present in 
real samples)^ and this increases the validity of the har- 
monic approximatioujiSSi Moreover, the infrared cutoff qc 
due to anharmonic effects depends on the applied strain 
(see Appendix |D]1 , decreasing as strain increases. This 
is consistent with sample to sample mobility differences 
of order unity recently reported in suspended monolayer 
graphene where strain u < 10"'* — 10""^ is naturally 
expected. 

Finally, we comment on a recent theory paper by Mar- 
iani and von Oppen^ where the T dependent resistivity 
of monolayer graphene has been fully discussed. In the 
high T regime, T 3> TbGj our results for the monolayer 
case agree with those of Ref. 0. However, at low T, i.e. 
T <C Tbg, the authors of Ref. found a new regime 
where the scalar potential contribution dominates. This 
happens because in Ref. the electron-phonon coupling 
from scalar potential is assumed to be much higher than 
the gauge potential coupling, unless T ^ Tqd, where 
Tod is the energy scale at which screening becomes rel- 
evant and scalar and gauge potentials become compara- 
ble. The new regime arises for ^ T ^ Tbg- With 
the parameter values used in the present work, however, 
the electron-phonon coupling due to scalar potential is 
always similar or smaller than the gauge potential (see 
Sec. nil Bp . Therefore, we can neglect this low T con- 
tribution since it gives higher power law behavior than 
the gauge potential. Recent T dependent resistivity due 
to in-plane phonons measured in single layer graphene at 
the high densities^^ seem to corroborate the latter pic- 
ture. 



VI. CONCLUSIONS 

In the present work we have studied the T dependent 
resistivity due to scattering by both acoustic in-plane 
phonons and FPs in doped, suspended bilayer graphene. 
We have found the bilayer membrane to follow the quali- 
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tative behavior of the monolayer cousini^ Different elec- 
tronic structure combine with different electron-phonon 
coupling to give the same parametric dependence in re- 
sistivity, and in particular the same T behavior. In paral- 
lel with the single layer, FPs dominate the phonon con- 
tribution to resistivity in the absence of strain, where 
a density independent mobility is obtained. This con- 
tribution is strongly suppressed by tension, similarly to 
monolayer graphenci^ However, an interesting quantita- 
tive difference with respect to suspended monolayer has 
been found. In the latter, as shown in Ref. |^, FPs limit 
room T mobility /i to values obtained for samples on sub- 
strate, ^ ^ 1 TD? / (Vs), when tension is absent. In bilayer, 
quantitative differences in electron-phonon coupling and 
elastic constants lead to a room T mobility enhanced by 
one order of ma gnitude, ^ 2Q mV(Vs), even m non- 
strained samples. This finding has obvious advantages for 
room T electronic applications. It has also been shown 
that for a correct description of acoustic phonon scat- 
tering in both monolayer and bilayer graphene, even at 
the qualitative level, coupling to both scalar and gauge 
potentials needs to be taken into account. 
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Appendix A: Collision integral 



The crucial step to get /kjg^^jj. is finding the scat- 
tering probability for a quasi-particle in state |k) to 
be scattered into state |k'). i.e. (since the pro- 

cess is quasi-elastic interband transitions are not al- 
lowed, meaning that both states belong to the same 
band). The scattering mechanism is encoded in the 
interaction 'Hint, which in the present case is given 
by Hep in Eq. (P^ . It is readily seen that scatter- 
ing occurs only through emission or absorption of one 
phonon or emission/absorption of two phonons. The 
initial and final states are thus tensorial products of 
the form \i) = |k) (X) |nq) or \i) = |k) (E) |nq,7iq/), and 
I/) = |k') ® |nq ± 1), I/) = |k') ® |nq ± 1, nq, ± 1) or 
I/) = |k') (g) |nq±l,n^, =Fl), where |7iq) and |nq,rtq/) 
represent one and two phonon states in the occupation 
number representation)^ and the electron like quasipar- 
ticle state is written according to the unitary transfor- 



mation in Eq. pi7| as |k) 



4lO)+e*H|0))/v^ 



(electron- hole symmetry guarantees the result is the same 
for both electron and hole doping). 



In order to obtain |(/|iJe 



with \i) and |/) as 



given above, we take the following steps, (i) Terms of the 
form V1V2, where Vi stands for scalar potential and V2 
for gauge potential induced matrix elements in Eq. (|23p , 
are neglected. It is easy to show that such terms come 
proportional to oscillatory factors e^'^^'' or e^'^^^' (in 
monolayer graphene e^'^'^ or e^"^^^' ) , stemming from the 



unitary transformation in Eq. (jlOp . These terms can 
safely be neglected in doing the summation over the di- 
rection of k and k' in the numerator of Eq. (|36|) . keep- 
ing 0k, k' fixed. The resistivity is then the sum of two 
independent contributions, originating from scalar and 
gauge potentials, well in the spirit of Matthiessen's em- 
pirical rulci^ (ii) For the scalar potential contribution 



the terms \ Vi^^e 



/4 are propor- 



tional to the overlap of states belonging to the same band 
(1 + cos26'k,k')/2) and can be written as 



The rate of change of /k due to scattering, the so- 
called collision integral /kj^^^jj appearing on the right 
hand side of the Boltzmann equation Eq. ([M|) . is the 
difference between the rate at which quasiparticles enter 
the state |k) and the rate at which they leave it. 



/k 



E 

k' 



[/k'(i-./k)w;^.-/k(i-/k')w^' 



(Al) 

where W/ is the scattering probability between state |i) 
and |/). Here we use Fermi's golden rule, which reads 



^ {f\mnt\i)^Si£f~£,), (A2) 



and is equivalent to rest upon the Born approximation for 
the differential scattering cross-section. In this appendix 
we provide explicit expressions for the integral collision 
Eq. (|Aip arising due to scattering by phonons in bilayer 
graphene (and also monolayer for comparison). 



v." 



1 -I- cos 26'k,k' 



(A3) 

The same manipulation holds for two phonon terms, with 
q, . (iii) For the gauge potential contribution 
there are oscillatory terms which, owing to the argument 
of point (i) , can be neglected. 



''2,q,k,k' 



-*(V+«k) 



2,-q,k,k' 



2 

— + — + /cA:'cose'k,k' ) ' (A4) 



where we used V to express the matrix elements given in 
Eq. for bilayer graphene without the term (TTk+TTk')- 
A similar manipulation holds for two phonon terms, with 

^2,q ^ '^2,q,q'- 

Finally, summing over phonon momenta and doing the 
thermal average, we can write as follows: when scat- 
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tering is via one phonon, 

= ^^Wi/(q,k, k')nq(5k',k+q'^(ek' - £k - ^^q) 
q 

+ ^ Wi/(q, k, k')(nq + l)(5k',k-q^(ek' - + fii^q), 

(A5) 

where the first term is due to absorption and the second 
to emission of a single phonon; when scattering involves 
two phonons, 

q.q' 

^k',k+q+q'(5(£k' - £k - flUJ^ - hhJ^,) 

q,q' 

^k',k-q-q'^(ek' - ek + ^l^q + ^l^q' ) 
27J- ^ 

+ — ^Wi.(q,q',k,k')(ri.q + l)nq/X 
q,q' 

^k',k-q+q'^(ek' - Ek + ?K^q - ?iWq,), (A6) 

where the first term is due to absorption of two FPs, 
the second to emission of two FPs, and the last one 
comes from absorption of a single FP and emission of an- 
other one. The kernels u;i,(q, k, k') and Wi?(q, q', k, k') 
represent the sum of the right hand side of Eq. (jA3[) 
with Eq. (jA4[) . as given in Eq. ([M)) . For Monolayer 
graphene take exactly the same formji only the ker- 
nels change, being given instead by Eq. ([55]) . The colli- 
sion integral may finally be put in the form, 

k' q,i^ 

X {[/k'(i - /kK - /k(i - /k')K + 1)] 

^k^k'+qf^lSk - £k' - fii^q) 
+ [/k'(l - /k)K + 1) - /k(l - /k')"q] 

^k.k'-qf^Cek - £k' + ^O) , (A7) 
for one phonon scattering processes, and 

A I scatt = J E E "'^^ ' ^ 

k' q.q' 

X {[/k'(l - ./k)(n-q + 1)(7V + 1) - /k(l - fk'hqriq' 

'5k,k'-q-q"5(ek - £k' + ^t^q + ^l^q' ) 
+ [./k'(l - /k)K + l)nq, - /k(l - /k')"q("q' + 1)] 

'5k,k'-q+q"5(£k - £k' + ^q - ^q' ) 
+ [./k'(l - /k)riqnq' - /k(l - /k')K + 1)K' + 1)] 

(^k.k'+q+q'^lEk - £k' - flUJ^ - hw^,) 
+ [/k'(l - /k)riq(nq. + 1) - /k(l - /k')K + l)"q'] 
^k,k'-Hq-q'<5(£k - £k' - fiWq -f- fiWq, ) } , (A8) 

for scattering through two FPs. 

Appendix B: Linearized collision integral 

In this appendix we derive the linearized version of the 
collision integrals given in Eqs. (jA7p and (|A8p . We start 



by expanding electron and phonon probability distribu- 
tions around their equilibrium values, 

/k = ' + <5/k, nq = + fcq, (Bl) 

where the variations can be written as (5/k = — gf^'fik 
[see Eq. ([26| ] and Snq = ~ g(yt^ ) Xq- The linearized col- 
lision integral J/k|g^^j.j is then obtained by expanding 
fkl up to first order in the variations i^^i'^^ 

I scatt ^ 

1. One phonon scattering 

This case follows closely the steps outlined in Ref. l39l . 
and for the case of monolayer graphene it has been de- 
rived in Ref. 0. Since the difference between monolayer 
and bilayer amounts to a different kernel ti;i,(q, k, k') in 
Eqs. (jA7|l . which does not play any role in the lineariza- 
tion, we can directly apply the result of Ref. to the 
present case. In order to set notation for the more elabo- 
rated case of two phonon scattering, we outline the main 
steps of the derivation in the following. 

We first note that at equilibrium detailed balance im- 
plies fl^^ I = 0, from which we get the relation 

/^^(l - fL"'W^^ = /f (1 - + 1), (B2) 

which can be easily verified by direct calculation)^ 
Therefore, in order to get the linearized collision integral 
it is enough to calculate the variation, 

S [/k'(l - /k)nq - /k(l - /k')K + 1)] = 

V t - Jk' ?^q + t -L - /k / 

(B3) 

The variations appearing on the right hand side of 
Eq. (|B3P can be computed easily by noting that 

^1-/^ (l-/(0))2 %n+iy (n(0) + l)2- 

Using Eq. (jB2p . and rewriting (5/k and Suq as 

^/k=/r(i-/r)^ ^nq=<n<^+i)^, 

''(B5) 

it is easy to show that 

S [/k'(l - /k)nq - /k(l - /k')K + 1)] = 
^/^(l - + 1) (V'k' + Xq - ^k) . (B6) 

The linearized collision integral may then be obtained, 
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^ E E "'-'(q' ^')fL'\^ - + 1) (^k' + Xq - V'k) <5k,k'+q<5(ek - ek' - K) 



hk^T 



k' q,i^ 



^ E E "'-(q- " /k°V^?^(<' + 1) (^k' - ^k - Xq) -Jk.k'-q^lsk - + K)- 

^ k' q,!/ 



r 



Now we introduce two typical approximations: con- 
sider phonons at equilibrium by taking Xq ~ 0, so that 



,(0) 



valid at not too low temperatures consider S [/k'(l — /k)(?^q + l)('^q' + 1) ^ /k(l — /k')'^q"q'] 



quasielastic scattering, with Ekj^k' 3> huq. The latter 
approximation enables us to rewrite the delta functions, 
'5(ek — Ek' ± ~^ S{ei^ — £k'): 8,nd to approximate 

/(°)(l-/(?))b5^ 



a j-(O) 



The linearized collision integral then reads 



2tt 



^^^.(q,k,k') 



k' q.i^ 



|^^(v:'k-^k')'5k,k'+q<5(ek-ekO, (B9) 



where we have used equalities w^{q,'k,'k') 



tOy(— q, k, k') and ujq 



-q- 



Finally, Eq. jM 



can be put in the form of Eq. ((27|). 

ALcatt=-E^'^^'^'(^''-^'''), (BIO) 
k' 

where T^k.k' is given in Eq. (|32p . 



/k 



/k 



1 - /k' 1 - /k "-q + 1 "q' + 1 



(B12) 



and 



'5 [/k'(l - ./k)(7^q + l)nq' - ./k(l - /k')"q("q' + 1)] 

(l-/f)(l-/^?^)(nr + l)(-i°^ + l)x 



"q ■ -/v-q' 

./k' nc^' /k ?iq 



1 - /k' ?^q' + 1 1 - /k ?lq + 1 



, (B13) 



the other two possibilities being related with these ones 
by a minus sign and k — k'. The variations appearing 
on the right hand side of Eqs. ()B12p and (jB13|) can be 
computed easily by using Eq. (jB4p . We then arrive at 
the variations 

S [/k'(l - ./k)(nq + 1)(?V + 1) - /k(l - /k')«q"q'] = 

(B14) 



and 



2. Two phonon scattering 

Now we proceed with the linearization of the collision 
integral in Eq. (jASp , originating from scattering processes 
involving two FPs. At equilibrium detailed balance is 
guaranteed, /k*'^ | s(.jjjj = 0, and the following two relations 
hold. 



AO) 



AO) 



.(0) 



AO) 

■Ik' 



AO) 

~ Jk' 

AO) 



AO) 

Jk 



(0) ^(0) 



+ 1 nq' + 1 ■ 



,(0) 



1 4-{0) (0) 



1 



(Bll) 



In order to get the linearized collision integral it is easy 
to see that we only need the following two variations. 



S [/k'(l - ./k)(nq + 1)?V - ./k(l - /k')«q("q' + 1)] 

a - + i)ny i^^^i±^^3iz:^ii^ 



ksT 



(B15) 



where we used Eq. (jBSP and the relations in Eq. (jBll[) . 

It is convenient to express the quantity (1 — /j^^') in 

terms of the difference (f^^ — fj^^ ) . For that we use the 
relation 



AO) AO) 

Jk' Jk 



1 - cxp[(ek' - ek)/fcsr] ' 



(B16) 



of Eq. (jB14p , where Ek' — ^k = ^q + ^^q' holds, we have 



which is easily verified by direct calculation. For the case 

e have 

(B17) 



1 



1 - exp[(ek' - eu)/kBT] 1 



15 



while in the case of Eq. (jB15|) , where et' ^ £k = fv^, 
fioj^i, we get 



nq{l + riq/) 



1 - exp[(ek' - ek)/fcsT] 



(B18) 



The variations in Eqs. (|B12p and (|B13p may then be cast 
in the form 



and 



S [/k'(l - /k)(?M + l)nq' - /k(l - /kO'^qK' + 1)] = 

p(o) .(o)x dn''q^ dn''°^ ^5?" (<<3k' + Xq' - <^k - Xq) 



(B20) 



S [/k'(l - /k)(nq + l)(nq- + 1) - /k(l - /k')"q"q'] = 



(0) 



.(0)^ 971^°^ 

"Jk / 



*^q°^ ksT (V?k' - <Pk - Xq - Xq' 



Inserting the latter results into the variation of the two 
phonon collision integral in Eq. (|A8[) . and recalling that 
there are two other terms which can be obtained from 
Eqs. (fBlQl) and (|B20| by multiplying a minus sign and 



(B19) 

substituting k — > k', we obtain 



scatt 



ci (0) a„(0) 

- -Tl^l^'^Fi'i^^'^^^nA' -/k 



d{njjjF) d[huj^,) l + nq + 



iq, 



((^Jk' - <y3k - Xq - Xq' 



k' q,q' 

'5k,k'-q-q"5(ek - Ek' + SWq + ^Wq, ) 

TTV^V^ , ,/N/.(0) .(0)x 5"-q^ ^"q°^ ^B?" . , n 

i"2^z^^f^(^''i''''^)(-^k' •^k )7i7i-^7i7i:^i — :— (</5k' + Xq' - v'k - Xq) x 



d{flhJ^) d{hhJ^,) Uq - nq/ 



k' q,q' 

<5k,k'-q+q'<^(ek - £k' + ftWq - f^q') 



k' q,q' 



(V'k' + Xq + Xq' - <y5k) X 



<5k,k'+q+q'^(£k - £k' - /^^q ^ ^q' ) 

k' q,q' q ^ ^ q / M m 

5\^M'+<i-ci'5{e]^ - ek' - + ftwq,). 



(B21) 



Now we introduce the two typical approximations: Wq = 
consider phonons to be in equilibrium, Xq « 0, so that with T'k.k' as given in Eq. ((33 
riq « n"^^ ] assume quasielastic scattering. From the lat- 
ter we obtain (5(ek — £k' ± ^q ± ^i^qO ^ '5(£k — Ck') 



It can be written in the form of Eq. (jBlOp , 



and /k' /k ~ -'-fe^ (^q ^q')' ^"^"^ ^'^^ linearized 
collision integral then reads, 



Appendix C: Calculating the resistivity 



6k 



271- rj^dhsr-i 

fcs^T— ((/3k' - Vk) X 



ft2 



5]«;j^(q,q',k,k') U 



F I F 
+ nq + Tlq' 



F F 



^|^'5k.k'+q+q'-5(.k-.k'), (B22) 



In this appendix we provide details regarding the 
calculation of the T— dependent resistivity for bilayer 
graphene. The variational method is used, with resistiv- 
ity given by Eq. ([55]) . Writing the numerator in Eq. ([55]) 
as in Eq. ([57|. the resistivity can be cast in the form 



Sn^e^kj, 



J dkdk' (K-u)^7'k,k' 



(CI) 



where we have used the fact that ?i;i?(q, q', k, k') is in- The remaining task is the calculation of the integral on 



variant under changes q — — q and q' 



-q', and the right hand side of Eq. (|C1 



16 



1. Scattering by in-plane phonons 



2. Scattering by flexural phonons 



This case follows closely the derivation to the Bloch 
law in 3-dimensional metalsi^ Inserting Eq. for 
Pk.k' into Eq. ([ClI) we get 



Vh 



j dkdk' (K • u)^ X 



^(ek-ek'), (C2) 



Inserting Eq. ^ for T'k.k' into Eq. (|CT]) we get 



^ WF(q,K - q,k, k') 



, driq, dfl°^ 

F , , ,F 



<5(ek - ek')x 



1 + nq + riK-q nq - riK-q 



, (C6) 



where we have already performed the sum over q. We can 
simplify the integral above by integrating over k and k' 
noting the presence of (5(ek~ek') and « —S{ep — e]^). 
The result reads, 



Vh / m \ 2 

J de^de^, (K • u)2 ^«;,(K, kpe^, fc^e^O^K^: 



where we have already performed the sum over q'. We 
can simplify the integral above by integrating over k and 
k' noting the presence of i5(ek — fk') and « —5{ep — 
Ek)- The result reads, 



WF (q, K - q, kpck, kpCk' ) -^-p — 

q q K-q 

/ F , F F F \ 

/ '^q + '^K-q '^q " '^K-q \ 



(C3) 



1 + "q + "K-q ?lq - 7lK-q 



(C7) 



with K — kp(ek — e^')- In order to proceed with the cal- 
culation we have to specify the kernel u'jy(K, fci?ek, kpe\^i) 
given in Eq. (|34p . Making use of the matrix elements in 
Eq. dSSl) we get. 



w^(K,fc_Fek,fc_Fek' 



[DUK/kF)\ 



2Vpvj 



The kernel z«F(q, K — q, fci?ek, fcFCk') is given by 
Eq. ([M]) with V"q V^^, (see Appendix O for a deriva- 
tion). Inserting the matrix elements in Eq. ([23]) it takes 
the explicit form 



Wir(q, K-q, kpek, kpe^') = wpiq, K, |K - q|) 
_ [D^{K/kF)YhS'k],\K-ci{ 



(C4) 



24V2mV24<^^-q| 



■, (C8) 



with Dg{x) as given in Eq. and where we have 

used the relation K ~ 2/;^ sin(6'k,k'/2)- The kernel de- 
pends only on 6'k,k', or equivalentlyX (the norm of K), 
as is the case of the rest of factors in the integrand of 
Eq. (1C3| but for (K • u)^. The latter can be written 
as (K • u)^ = cos^ 7, and the angular integration is 
then conveniently done by integrating over 7 keeping 
6'k.k' = 6v — = 6 constant, and integrate over after- 
ward, or equivalently K. Using dO = dK/^^kp — /\^/4, 
the resistivity becomes 



2h^e^k 



2fcf 



dK 



(0) 
K 



(C5) 



Inserting Eq. (jC4p for the kernel {K) into Eq. (jCSP we 
readily obtain Eq. ([55]). From that it is straightforward 
to calculate analytically the two limiting cases, T <C Tbg 
and T^Tbg-^ 



with (x) as given in Eq. (j50p , and where we have used 
the relation K ~ 2kp sin(6'k,k'/2) and assumed Wq given 
by Eq. ([5]). In deriving Eq. (|C8|) we used cos'^{4> — 4>') = 
[1 H- cos(2(/) — 2(/)')]/2 and dropped the oscillatory part. 
The sum over q can be replaced by an integral, — > 

J qdqdcj), and owing to the relation = |K — qp = 
+ — 2qK cos (f> we can write the resistivity as 



K+q\ 



dQ 



de^dOi,, (K • u 
QivF{q,K,Q) 



,2 r, dn, 

d^n 



dnq 



\K~-q\ 



q^K^ - + g2 _ g2 /4 Q 



where 



we 



1 + + ng 
used 



F F 
^q 



, (C9) 



dQQ/dq^K^- {K^ + q^ -Q^) I A. As in Sec. [CTl the 
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angular integration over 0k and 0k' is conveniently done 
by integrating over 7, with (K • u)^ = cos^ 7, keeping 
^k.k' = ^'k — ^k' ^ and (7 and |K — q| = Q constant, 
and integrate over afterward, q and Q. The resistivity 
may then be written as 



\K+q 



dq 



q^ duq 



dQ- 



Vk'p - Kyi 



K-q\ 



from which we obtain 



2K 



dqqnq{nq+l) — 



The resistivity may finally be written as a simple integral 
over K, 



■In 



(ML 

\ haq^ 



q^K^ - K 




'''dKSmk^K%,in, + l), (C12) 







kj, - A'2/4 



form which Eq. (|49p is readily obtained. 



where d9 = dK/ \/kp — if 2/4 has been used, and we used 
Eq. ([CS]) for the kernel. 



Non-strained flexural phonons 



In the absence of strain the FP dispersion reads 
ujq = aq^ . After rescaling momentum as a; — > i = 
x{ha/kBTY^^ we can rewrite the resistivity as, 



{ksTy 



2^TT^he^p^vj,kj,a'^ 



2kf 



f 

Jo 



dqqnq{nq+l) 



K+q\ 



dQ 



\K-q\ 



dk'mm^. 

kj,-K^/4 



Strained flexural phonons 



The flexural phonon dispersion in the isotropic approx- 
imation is UJ^ sa \Ja^q^ + uv\q^ [see Eq. ([8])]. After 

rescaling momenta a; — > i = xhvLU^^'^ I {ksT) the re- 
sistivity in Eq. (jClOP takes the form given in Eq. ([55)1 . 
The low T regime is detailed in the main text. Here we 
concentrate in the high T regime, showing in particular 
how to obtain Eq. ([60]) for the integrals over q and Q in 
Eq. dSH). 

We start by writing the Q integral in Eq. ([55]) as 



Al,K,q) 



1 



\K+q\ 



dQ 



\K-q\ 



q^k^- [k^ + q^-QA /4 



1 + llq + rig Uq — Uq 



(Cll) 



7'Q^ + Q2 



^^2qi + g2 + ^^2Q4 ^ Q 
I + Uq + Uf^ 



7'Q^ + 



The integral over Q is infrared divergent, and is thus 
dominated by the contribution K ^ q. Defining the small 
quantity 5x = \K — q\, and noting that for Q 1 we have 
Uq ^ 1/Q^ ^ 1, it is possible to identify the dominant 
contribution in the Q integral as. 



\K+q\ 



dQ 



K-q\ 



Q' 



q2j^2 _ 

~2 r\2 \ 

q -Q \ 



1 



2K' 



2K n - + 
dQ^ 



2K 

5x 



It is now obvious that the q integral has a logarith- 
mic divergence for q ^ K. Note, however, that in the 
present theory phonons have an infrared cutoff, so that 
min jif — g| = where q^ ^1 \s either due to strain or 
anharmonic effects. The dominant contribution to the q 
integral is then coming from the maximum of l/j/f — g|. 



(C13) 

with 7 = \/2L0q^ /ujq, Haviug in mind that high T im- 
plies K <^ 1, we consider the integration in Eq. (jC13[) in 
two limiting cases: when q ^ k <^ \ and for q ^ k. 
In the former case, since 5 <C 1 and Q ^ 1 hold, we 
can linearize the dispersion relation and approximate the 
Bose-Einstein distribution function by Uq sa 1/q and 
rig « 1/Q (as discussed in the main text, finite strain 
implies kp ^ <?*, so that the linearization of the fiexu- 
ral phonon dispersion can be taken when q < K <^ I 
holds). The integral over Q in Eq. (|C13P may then be 
approximated by 



2^(7, K, q) 



K+q\ 



dQ 



2qQ 



\K~q\ 



q^K^ -[K^ + f-Q'^) /4 

(C14) 
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and the integral can be done as 

4gQ 



\K+q\ 



clQ 



K-q\ 



2q arctan 



where we have defined 



\K+q\ 



= 27rq = I{q), 



\K~q\ 



(C15) 



Y{q,Q,K) = -{q-K-Q){q-K+Q){q+K-Q){q+K+Q). 

On the other hand, for q ^ K the integration region is 
concentrated around q. We may then write the integral 
in Eq. (jC13P as a slowly varying function, which we can 
take out of the integral, multiplied by an integral of the 
form of that in Eq. (|C14p . 



I{l,K,q) 

\K-q\ 



1 + 2nq 



nq{nq + 1) 



dQ 



Q 



Iq^K^- i^lO + q^^Q^) /4 
TTq^n.ging + 1) ( 2/7^*T^ , 1 



1 + 2nq nq{nq + 1) 



(C16) 



Since for q < K the later result reduces to 27rg, as in 
Eq. (jClSp . we can use 1(7, q) in Eq. (jC16P to approximate 
the Q integral, Eq. ([C13l) . in the full region g < A' <C 1 
to g ^ K. This has been tested numerically to be a 
good approximation as long as X ^ 1. The q integral 
in Eq. (|56p may then be cast in the K independent form 
given in Eq. ([50]) . 



Appendix D: Perturbative treatment of anharmonic 
effects 



Recently it has been demonstrated how anharmonic 
effects in stiff membranes as graphene are highly sup- 
pressed by applying tension)^ In this Appendix we show 
how the infrared cutoff q^ of our (harmonic) theory de- 
pends on the applied strain. This is consistent with the 
sample to sample differences of order unity reported in 
recent measurements of electron mobilities in doped sus- 
pended monolayer graphenci^ As argued in Ref. |6|, com- 
paring qc ~ O.lA 1^ (the infrared cutoff for the har- 
monic theory in the absence of strain), and q* [the mo- 
mentum scale associated to the presence of strain, de- 
fined by Eq. ([55]) ] gives u ~ 10"'' — 10"'^ as the strain 



involved in such kind of experiments. In order to esti- 
mate properly this number we have to know how q^. is 



+ 



Figure 3: First order diagrams for E (q). The curly line rep- 
resents the free correlator of flexural phonons, the dashed 
line corresponds to the 4-leg interaction vertex associated to 
Rij,ki (k). The second diagram is 0, since the q = Fourier 
component of the transverse projector is integrated out during 
the in-plane modes Gaussian integration. 



affected by strain. The discussion can be applied to both 
cases, monolayer or bilayer graphene. 

In order to estimate anharmonic effects in the disper- 
sion relation of FPs, since Eq. ([1]) is quadratic in the 
in-plane displacements we can integrate them out to ob- 
tain the effective free-energy for the out-of-plane degree 
3m^ 



•If 



1 



dxdy (y h) 



dxdyR,j^kid'hd^hd^hd^h, 
(Dl) 



where the four-point-coupling fourth-order tensor can be 
written as Rau = ^PijPki, the operator is the 



Hj,kl 

transverse projector = (^^) ^ £ik£jid''d\ and Kq = 
4t4p+\) ^ order to include the effect of strain, we add 
to Eq. (jPip the simplest term which breaks rotational 
symmetry 

^7 / dxdy{\Ihf , 



where 7 is a sample-dependent coefficient with units of 
tension which can be related with the strain of the sam- 
ple. This approach follows the spirit of the effective 
isotropic dispersion relation of FPs introduced in Eq. ([5]) . 
If we add to Eq. ([T]) the most general term to first or- 
der in the derivatives of the displacement fields which 
breaks rotational symmetry and then we integrate out 
the in-plane degrees of freedom, we obtain a new two- 
point vertex whose contribution to the renormalization of 
the bending rigidity is weak and can be neglectedi^. We 
are going to study the Fourier component of the height- 
height correlation function 



G(q) = (|Mq)l 



m(q)|/i(q)|2e^^"[''('l)l. 



where obviously Z = j Vh (q) e"^ ^' is the partition func- 
tion of the system, the effective action is nothing but 
S''f^[h{qj\ = -liF'^ff [h{qj\, and the Fourier trans- 
formed effective free energy reads 
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T'i^ [/i(q)] = 



(27r)^ 



4 , 7 2 

— (7 H — q 
2 2 



IMq) 



i2 

n 7^^91929394%.^; (qi + q2) -^^'^ (qi + q2 + qs + q4) , (D2) 
^=1 (27r) 



r 



where R^^u (k) 



transformed 

„k 1,1 



^Pf^ (k) i^^ (k), and the Fourier 
reads 



transverse projector 



k~'^eikejik^h} = — -^^^ It is important to note that 
the q = Fourier component of the transverse projector 
is integrated out during the Gaussian integration of the 
in-plane modes4^ In the harmonic approximation and in 
the absence of strain (7 = 0), the (free) correlator is given 

by 



g(0) 



KbT 
Kg4 



(D3) 



When we assume a quadratic dispersion relation we are 
taking this correlator as the proper one. This approx- 
imation is obviously affected by the presence of strain 
and anharmonic effects (also affected by strain), which 
renormalizes the bending rigidity k. Then, we can write 
(q) cx K (q) q** and study the renormalization of k 
from the Dyson equation 



(q) = (g(°) (q))"'+S(q)^ 



(D4) 



where now the correlator in the harmonic theory, includ- 
ing the effect of strain, is given by 



G(°) (q) 



KbT 



(D5) 



In order to estimate the anharmonic effects we compute 
the first order diagrams for self-energy, showed in Fig. [31 
Only the first diagram gives a non-zero contribution 



Follow the Ginzburg criterion;^ we estimate the cutoff of 
the theory nothing but comparing each correcting term 
of Eq. (jD7P with the bare value of k. As we have already 
mentioned, there are two different cutoffs, the one given 
by strain in the harmonic approximation, and the other 
one associated to anharmonic effects. The first one is 
given by the second term of Eq. (|D7P , and it is nothing 
but q* . 



(D8) 




2.10-" 



4.10"" 



10" 



Figure 4: The infrared cutoff as a function of the ap- 
plied strain u for monolayer (in red dashed Une) and bilayer 
graphene (in blue). In both cases T — 300 K. 



(q) = 4/3 



(27r) 



2 9 r 9 q ri ij 



R,,,ki (k)G(") (q-k). 



Replacing this result in equation Eq. (|D4|) we obtain 

2 



(D6) 



K (q) = n - 



7 



KbTKq 



1 (q,:'^) 



(27r)^ K|q-k|4 + 7|q-k|2- 
(D7) 



Identifying this result with the momentum scale defined 
by Eq. (j55p . we deduce the relation between 7 and u: 
7 = (A 2^i)u. 

The cutoff of the harmonic theory is also affected by 
strain. Following the same criterion, its value comes from 
the solution to the equation 



{2ttY k Jo Jo (k^ + q'^ - 2kq cos (6*))^ + {q*y (fc2 + q^ - 2kq cos (9)) 

I 



(D9) 



In the absence of strain, we have q^ = J ^^iJJ" , which gives at T = 300 K the value q^ = 0.178A ^ in the case 
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of monolayer, and Qc = 0.126A in the case of bilayer creases, so the unavoidable little strain present in real 
graphene. Its dependence on the applied strain is shown samples increases the validity of the harmonic approxi- 
in Fig. 01 It is clear that Qc decreases as the strain in- mation. 
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